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摘 要 为 了 增加 间断 Galerkin(Discontinuous Galerkin, DG) 方法 在 非 定常 流动 中 的 求解 效率 , 本文 开展 了 非 定常 流 
动 的 隐 式 DG 方法 研究 . 隐 式 DG 方法 的 构造 采用 二 阶 向 后 差分 格式 (BDF2) 进行 时 间 项 离散 ， 非 线性 代数 系统 的 求解 
基于 Newton 和 迭代 法 , 采用 块 对 称 Gauss-Seidel(SGS) 迭代 法 对 线性 方程 组 进行 了 求解 . 基于 所 发 展 的 非 定常 流动 的 隐 式 
DG 方法 , 分 别 对 等 粹 圆柱 扰 流 和 卡门 涡 街 (Re=100) 现象 进行 了 数值 模拟 . 研究 结果 表明 , 所 发 展 的 隐 式 DG 方法 能 够 
达到 设计 精度 ， 能 够 在 高 出 显 式 方 法 两 个 数量 级 的 时 间 步 长 上 保持 稳定 ， 具 有 高 的 求解 效率 ， 且 计算 结果 与 显 式 方法 和 相 
关 文 献 均 吻 合 较 好 . 
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0 引 As 

高 精度 CFD 方法 是 空间 高 于 二 阶 精度 扒 方 程 
离散 方法 ， 相 比 于 低 阶 方法 ， 由 于 具有 在 较 低 的 计 
算 成 本 上 得 到 更 高 计算 精度 的 潜在 特性 而 受到 了 关 
注 趾 , 目前 发 展 了 一 系列 的 高 精度 CFD 方法 , 例如 
高 阶 有 限 体积 、DG 方法 、 谱 方法 、 谱 差分 / 谱 体积 、 
通 量 重 构 等 . 其 中 DG 方法 具有 高 阶 精 度 、 易 于 构 
造 、 播 值 模板 小 、 灵 活 处 理 间 断 问题 和 易于 实现 并 
行 计算 等 特点 ， 自 从 上 世纪 80 年 代 未 Shu 等 剖 提 
出 DG 方法 求解 守恒 律 问题 的 基本 理论 以 来 便 得 到 
了 大 量 研究 8-5. 

DG 方法 结合 了 有 限 元 方法 (FEM) 和 有 限 体积 
方法 (FVM) 的 特点 . 在 单元 采用 与 FEM 相似 的 多 
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项 式 允 近 集 略 ,通过 增加 基 活 数 个 数 以 提高 计算 精 
度 ; 不 同 点 在 于 其 在 每 个 单元 上 的 多 项 式 分 布 是 完 
全 独立 的 ， 且 网 格 单元 界面 上 人 允许 间断 产生 . 单元 
之 间 的 信息 传递 采用 了 FVM 中 Riemann 求解 属 计 
算 界面 通 量 的 方法 ,保持 了 格式 的 迎风 特性 ， 且 在 
超声 速 流 场 中 限制 器 的 构造 完全 继承 于 FVM. 
近年 来 , 国内 外 学 者 针对 于 DG 方法 的 构造 进 
行 了 大 量 的 研究 . 张 来 平 日 对 于 高 精度 DG/FV 混 
合算 法 进行 了 研究 ， 其 在 DG 框架 下 采用 重 构 思想 
通过 扩大 插值 模板 在 存储 其 基本 不 变 的 前 提 下 到 达 
更 高 阶 精度 。 于 剑 、Jiang 等 "3 基于 四 边 形 网 格 
分 别 研究 了 层 流 流动 DG 方法 的 构造 和 非 定常 流动 
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的 隐 式 求解 方法 . 部 海兵 加 基于 Gauss-Seidel 迭代 
研究 了 无 务 流 动 的 隐 式 DG 方法 。Wangt 详细 阅 
述 了 非 定常 无 黏 流动 的 隐 式 DG 方法 及 p 型 多 重 网 
格 方法 的 构造 。Perssonl11 针对 隐 式 DG 方法 中 病 
态 线性 方程 组 的 预 处 理 方法 和 迭代 方法 进行 了 详细 
研究 . Patrickbd 采用 解析 方法 推导 了 无 黏 流 动 隐 式 
DG 方法 Jacobi 矩阵 的 计算 . 谭 勤 学 上 3] 对 雷诺 时 均 
N.S 方程 (RANS) 采用 DG 方法 离散 , 通过 数值 模拟 
研究 了 燃气 轮机 空气 系统 中 的 流动 和 传 热 问 题 ， 研 
究 表明 DG 方法 的 结果 更 接近 于 实验 值 . 此 外 , DG 
方法 在 压气 机 等 部 件 的 内 部 流动 研究 中 也 得 到 了 相 
应 的 应 用 1s,19 

研究 者 们 采用 DG 方法 进行 流 场 计算 时 广泛 采 
用 了 显 式 Runge-Kutta 方法 对 时 间 项 进行 离散 ， 通 
过 时 间 推 进 方 法 得 到 收敛 的 流 场 结果 。 虽然 此 方法 
简单 旦 易于 构造 ,但 是 由 于 数值 稳定 性 严重 限制 了 
时 间 步 长 使 得 流 场 求解 效率 异常 低下 . 例如 , 采用 2 
阶 精度 的 DG 方法 进行 流 场 求解 时 保证 绝 大 部 分 算 


其 中 g = [wu; slT， wi; 和 g; 分 别 表 示 第 ;个 基 函 数 和 
其 系数 ，N, 为 每 个 单元 上 的 自由 度 总 数 . 基 函 数 的 
选取 可 采用 多 种 形式 , 例如 Lagrange 基 薄 数 和 任何 
一 种 标准 正 交 基 函 数 . 

其 次 , 将 式 (3) 代入 式 (1), (2) 中 , 基于 加 权 余 
量 法 ( 变 分 原理 ) 使 余 量 Rh 正 交 于 基 沁 数 空间 ， 即 
基 函 数 与 余 量 方程 在 每 个 控制 体 上 的 积分 为 零 : 

i: R= Ty Fin (Wn)—V. Fvis(up, Vup) 
(4) 
| pi RndV =0 1l1<i< Np 
DE 


将 其 展开 便 得 到 了 DG 方法 的 半 离 散 方 程 : 


1 4 多 = VD0d5 -/ UVPpidf 


piundo= 上 (Finv — Fvis) . Vpidf2— (5) 
DE D 


k 


全 区 的 CrXn015, 在 更 高 阶 本 下 这 一 限 抽 人、 多 wmv ~ Pus) mrds 


更 加 严 奇 19. 因此 提升 流 场 求解 效率 是 DG 考区 
行 工程 应 用 研究 亟待 解决 的 关键 问题 之 

在 流动 仿真 中 提升 求解 效率 的 3 有 隐 式 
计算 方法 、 并 行 计算 和 多 重 网 格 徐 % 践 中 隐 式 方法 
可 论 上 不 妥 时 间 步 长 的 限制 (和 亲 件 和 定性， 
可 显著 提升 流 场 求解 效率 . 丹 此 ' 本 文 开展 非 定常 层 
流 流动 的 隐 式 DG 方法 研究 , 基于 所 发 展 的 隐 式 DG 
方法 , 分 别 对 等 炉 贺 柱 扰 流 和 卡门 涡 街 (Re = 100) 
现象 进行 了 数值 模拟 ,以 验证 本 文 所 发 展 的 非 定 党 
层 流 流动 的 隐 式 DG 方法 的 求解 精度 及 稳定 性 . 


1 DG 方法 与 FVM 对 比 


1.1 间断 Galerkin 方法 AS 
非 定常 可 压缩 N-S 方程 为 : BS 
Ou 
可 +V:. Finv(u)= Vv. Fvis(u,Vu) (1) 


其 中 ， 人 为 守恒 变量 ， Finv 和 下 vis 分 别 为 黏 性 通 量 
项 和 无 黏 通 量 项 。 其 中 黏 性 通 量 项 的 本 构 关 系 满 足 
Stokes 假设 , 热流 量 满足 Fourier 定律 . 为 了 采用 DG 
方法 进行 方程 离散 ,引入 辅助 变量 s: 
s=Vu (2) 
首先 ,将 求解 区 域 进 行 单元 剖 分 ， 在 每 个 单元 
上 将 守恒 变量 v 和 辅助 变量 s 采用 多 项 式 通 近 : 


Np 
q = qi(t)pi(2) (3) 


其 中 w* 和 [(Fiwy 一 下 vis). n]* 为 数值 通 量 , n 为 界 
面 外 法 向 向 量 。 本 文中 无 莫 项 通 量 和 香 性 项 通 量 分 
别 采 用 Lax-Friedrich 格式 和 BR1 格式 。 
最 后 , 将 式 (5) 时 间 项 和 空间 项 分 离 , 简化 得 到 
< d 
u 


2 My = R(u) (6) 


1.20 和 潮 体 积 方法 


\ FVM 方法 的 解 在 单元 上 采用 常数 近似 , 根据 数 


据 存储 点 与 网 格 节 点 位 置 关系 通常 有 单元 顶点 法 和 
单元 中 心 法 两 种 . 控制 体 中 解 的 分 布 为 : 


1 
7k 


式 中 V* 为 控制 体 D* 的 体积 (面积 ). 
同样 , 首先 将 其 带 入 余 量 方程 (4) 中 的 i 式 , 并 
在 单元 上 积分 : 
十 V (Finv(U) — Fvis(U, Vu)) dV =0 (8) 
ps Ot 
其 次 , 利用 Gauss 定理 将 散 度 的 体 ( 面 ) 积分 变 
为 界面 上 的 面 ( 线 ) 积分 : 


dx 1 


a (F(T)— Fi (U, Va)): nds=0 (9) 


V1S 


4 期 


式 (9) 则 为 FVM 的 最 终 离散 格式 。 提 高 空间 
逼近 精度 的 关键 在 于 重 构 守 恒 变 量 在 单元 上 的 分 布 ， 
格式 的 精度 将 由 重 构 界 面 通 量 的 精度 决定 。 

相 比 于 DG 方法 的 离散 过 程 , FVM 在 离散 时 相 
当 于 DG 方法 中 守恒 变量 yp = 1 时 的 情形 . 因此 , 在 
离散 方法 上 FVM 可 以 认为 是 DG 方法 在 基 函 数 为 
1 时 的 特例 . 


2 非 定 常 隐 式 DG 方法 构造 


显 式 DG 方法 在 计算 式 (6) 右 端 残 量 项 时 采用 
当前 已 知 时 间 层 上 的 守恒 变量 值 R(w"). 而 隐 式 DG 
方法 是 指 残 量 项 基于 待 求 时 间 层 n 二 1 时 刻 的 守恒 
变量 值 : 


(w"™) (10) 


显 式 方法 和 隐 式 方法 的 本 质 差异 在 于 右 端 残 量 项 的 
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首先 ，Jacobi 矩阵 dR/dup 的 计算 对 于 条 性 流 

动 若 采用 解析 方法 通常 是 困难 的 ， 因 此 本 文采 用 数 

值 方法 进行 近似 计算 : 
dR 


R(ur +e)— 
了 |v =ur 
u 


€ 


R(ur) (15) 
其 中 s 为 一 小 量 , 例如 可 采用 eX ||ul| x 10-s。 

其 次 ,对 于 迭代 过 程 (14) 中 并 式 的 求解 . 对 其 
进行 整理 可 写 为 : 


A Cs 加 证 | =B (16) 
3M dR 
2At dan ” 
M /3 ,01 1 
B 一 一 A (Be — 2up 十 3 uh 
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三 ”计算 是 基于 哪 一 时 间 层 上 守恒 变量 值 的 ， 两 者 与 时 < Te) 
LO 间 项 离散 格式 无 关 ， 因 此 任何 一 人 SS 
G) ”离散 格式 均 可 用 来 构造 隐 式 方法 。 本 文 将 以 时 间 3 (16) 为 大 型 线性 方程 组 , 待 求 值 为 网 格 单元 
中 。 二 阶 向 后 差分 格式 (BDF2) st Us 由 度 上 守恒 变量 的 值 .然而 其 系数 矩阵 A 
@ 方法 构建 , 待 求 方程 为 ， 8 于 通常 具有 很 大 的 条 件数 , 造成 了 方程 (16) 严重 病态 ， 
LO M (Ban ut 1 滞 - 20 (11) 这 为 方程 组 的 直接 求解 带 来 了 困难 . 因此 , 在 CPD 
忆 At \2 I 中 的 隐 式 方法 广泛 采用 对 系数 矩阵 先 预 处 理 再 进行 
非 定常 隐 式 方法 构建 的 3 2 (11) 精确 os 国人 Gauss-Seidual 
| 。 成立 , 其 所 有 的 工作 均 是 帮 儿 各 计算 效率 /存储 量 同 。 顶 处 理 、LU 分 解 党 ”。 本 文采 用 块 Gauss-Seidual 
:"” 式 (11) 精确 成 立 之 间 的 矛盾 . 由 于 式 (11) 经 离散 后 ss 程 (6 的 求解 的 , 求解 过 程 均 是 以 每 
:三 。 ”为 一 大 型 非 线性 代数 系统 ， 为 了 对 其 进行 求解 本 文 Me 度 的 守恒 变 基 值 构成 的 矩阵 为 元 
全 。 采用 Newton 迭代 法 ， 人 Us 
个 wntt 到 M 3 nt1 有 也 十 到 D+ U 
二 Wy ( ee 和 D)Awn+15=B-U.Auntl-l j=1,3,5. 机 
© —R(un'')=0 内 CD+D)Aw 一 一 开 Au j=2,4,6. 
人 式 中 ， 厂 D, U 分 别 为 下 三 角 阵 、 对 角 阵 和 上 三 角 
dR : > n+1,7 Wr tl wnt 二 
Te (wnt ynthil) ~ R(t .> (13) 。 阵 . Av = 一 他 。 在 进行 第 奇 次 和 
0 % 第 偶 次 远 代 时 分 别 采用 向 下 扫 措 和 向 上 打 措 依次 求 
式 中 j = 1,2,… ,6, 6 为 克 代 总 次 数 ，w?117 为 第 解 ; 每 次 只 需要 对 每 个 单元 的 系数 矩阵 进行 求 逆 运 


j 次 迭代 的 待 求 时 间 层 n+1 上 的 守恒 变量 值 . 将 式 
(12) 带 入 式 (13) 中 ,可 得 到 迭代 过 程 为 : 


i: wnt = wu? 

sie 3 M dR nn 二 1,7 n+1,7—1 

i 了 Wh ) 
M 3 n+1,7—1 n 生 二 


Rar)) 


ntl wnt 


11i: Up 


算 , 显著 降低 了 条 件数 和 计算 量 . 

至 此 便 完成 了 对 非 线 性 代数 系统 式 (11) 的 求解 
过 程 。 对 于 非 定常 流动 ， 时 间 推 进 的 每 一 个 时 间 步 
上 式 (11) 的 计算 准确 度 均 与 流 场 精度 密切 相关 ; 另 
一 方面 , 由 于 对 式 (11) 的 计算 采用 了 达 代 求解 策略 ， 
总 是 存在 误差 使 得 其 不 可 能 精确 成 立 。 因 此 ， 对 于 
非 定常 流动 的 适 代 收敛 判 据 是 在 保证 较 少 的 迭代 次 
数 与 式 (11) 计算 准确 度 之 间 做 出 的 折 中 : 


4 1 也 一 nn 
| 芝 (Be pt 路- Rlunt!)l<e lull (18) 
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本 文中 上 = 10-4 9, 圆 为 远 场 边界 条 件 、 内 部 圆 为 无 穿 透 (反射 ) 边界 


需要 指出 的 是 : 1) 非 定常 流动 中 时 间 步 长 At 必 
须 在 所 有 网 格 上 保持 一 致 ; 2) 式 (18) 为 隐 式 非 定常 
流动 判断 Newton 和 迭代 收敛 的 唯一 判 据 ， 因 此 在 计 
算 过 程 中 可 以 采用 与 隐 式 定常 流动 中 一 致 的 加 速 方 
案 , 例如 Jacopi 矩阵 延迟 更 新 、 线 性 方程 组 (16) 有 
限 步 从 代 等 . 3) 隐 式 DG 方法 虽然 能 够 显著 提升 时 
间 步 长 ， 但 是 在 每 个 时 间 步 上 通常 需要 数 倍 于 显 式 
方法 的 计算 量 与 存储 量 ， 尤 其 是 在 高 阶 精度 和 大 网 
格 量 的 情形 下 。 所 以 ， 想 要 彻底 改变 隐 式 方法 的 计 
算 效 率 寻 求 更 好 的 计算 策略 至 关 重要 . 

3 数值 结 

基于 本 文 构造 的 非 定常 隐 式 DG 方法 的 计算 方 
案 ， 对 无 务 圆 柱 扰 流 和 卡门 涡 街 现象 (Re=40) 进行 
了 数值 模拟 研究 . 前 者 由 在 证 明 本 文 计算 方案 的 正 
确 性 ,后 者 通过 非 定常 的 卡门 涡 街 现象 以 验证 非 定 
常 隐 式 DG 方法 的 计算 效率 和 计算 精度 . 
3.1 无 黏 圆柱 扰 流 


对 称 的 特点 使 得 流 场 结构 理论 上 高 度 对 
的 壁面 利于 研究 高 精度 CFD 方法 天 度 的 依 


赖 性 . 因此 在 网 格 量 分 别 为 74、59& 品 3000 的 三 套 网 
情 上 进行 了 2 阶 精度 (1) < 度 (ps) 的 无 攻 
圆柱 扰 流 流 场 计算 . 这 


图 1 网 格 示意 图 (74, ps) 
Fig. 1 Grids of isentropic How around cylinder (74, p3) 


圆柱 扰 流 来 流 条 件 为 Ma = 0.38, pw = 101325 
Pa 7 = 273 K; 高 阶 壁面 的 构造 基于 圆 弧 通 近 ; 
网 格 量 为 74 的 计算 域 示 意图 如 图 1 所 示 ; 外 部 


2 
夸 昌 


条 件 。 图 2 为 网 格 量 为 3000 的 网 格 上 2 阶 精度 
(Pp1) 时 显 式 方法 与 非 定常 隐 式 方法 的 残 差 收敛 历史 
图 ， 如 图 所 示 ， 隐 式 方 法 的 计算 效率 显著 高 于 显 式 
方法 . 图 3 表示 圆柱 表面 压力 系数 分 布 , 如 图 所 示 ， 
圆柱 上 下 壁面 高 度 重 合 且 左右 对 称 分 布 符合 无 秋 贺 
柱 扰 流 特点 ， 另 外 其 与 Bassin9 的 计算 结果 高 度 吻 
合 也 证 明了 计算 的 准确 性 . 图 4 和 图 5 分 别 表 示 
4 阶 精度 时 马赫 数 Ma 和 密度 分 布 等 值 线 图 ， 其 值 
分 别 为 Mai = 0.036 x (i 一 1),i = 了,26;pi = 
0.02 x (i—1)+0.94,i= 1,2,... ,23, 


一 -一 JImplicit (p1) 
一 一 一 Explicit (pi, At=10-9) 


ImplicitDG (p;) 
F.Bassi'sresult (p;) 


图 3 上 、 下 表面 压力 系数 分 布 图 


Fig. 3 Pressure coefficient distribution around cylinder 


表 1 为 2 阶 精度 至 4 阶 精度 时 不 同 网 格 量 上 得 
到 的 炉 误差 . 如 表 所 示 , 在 3 阶 和 4 阶 精 度 时 流 场 
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的 真实 精度 均 达 到 了 设计 值 ; 但 2 阶 精度 时 真实 计 
算 精度 低 于 设计 精度 ,这 是 由 于 在 2 阶 精度 时 单元 
的 等 参 变换 与 仿 射 变换 相同 ， 使 得 边界 单元 采用 了 
直 边 界 而 造成 的 . 这 表明 了 高 阶 壁面 近似 对 高 精度 
CFD 方法 至 关 重 要 ; 另外 , 本 文 构建 的 隐 式 DG 方 
法 实际 精度 能 够 达到 设计 精度 ， 


图 4 Ma 分 布 云图 (3000, ps ) 
Fig. 4 Ma contour around cylinder (3000, < 


xX~ 


图 5 密度 分 布 云 图 (3000, p3) 
Fig. 5 Density contour around cylinder (3000, p3) 


表 1 基于 非 定常 隐 式 DG 方法 的 精度 验证 
‘Table 1 Accuracy verification using implicit DG 


Element D1 p2 D3 


number Error Order Error Order Order 
74 1.4x10-? —- 4.1x10-3 — 1.1x10-3 一 
528 3.8x1073 1.33 1.0x10-4 3.75 3.2x10-6 5.87 


3000 7.1x10-4 1.49 3.1x10-6 4.06 1.0x10-8 5.90 
3.2 卡门 涡 街 

卡门 涡 街 现象 数值 模拟 的 计算 域 如 图 6 所 示 ， 

其 网 格 量 为 6368, 圆柱 壁面 分 布 40 个 网 格 点 , 同样 


Error 


G \\ Time step(implicit) 
四 
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地 对 其 边界 单元 进行 数值 积分 时 采用 圆 弧 允 近 边界 
边 的 曲 边 三 角形 进行 计算 。 自 由 来 流 参 数 为 pe = 
101325 Pa, Ts = 273 K, Ma = 0.3， 以 圆柱 直径 和 
远 场 边界 参数 定义 的 雷诺 数 Re=100; 圆柱 壁面 为 等 
温 (Twan = 273 K) 无 滑 移 边 界 条 件 。 本 文 对 其 进行 
了 2 阶 精度 (p1) 至 4 阶 精 度 (ps) 的 计算 . 所 有 精 
度 下 的 流 场 初始 化 参数 均 为 同一 个 收敛 的 2 阶 定常 
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图 6 网 格 示意 图 


RS Fig. 6 Grids of Karmann vortex street phenomenon 


表 2 为 不 同 精度 下 隐 式 DG 方法 的 时 间 步 长 与 
经 测试 得 到 保持 稳定 的 显 式 方法 时 间 步 长 对 比 。 如 
表 所 示 ， 隐 式 DG 方法 的 时 间 步 长 通常 比 显 式 方法 
高 出 两 个 数量 级 , 这 表明 本 文 发 展 的 隐 式 DG 方法 
能 够 在 远大 于 显 式 方法 的 时 间 步 长 下 保持 稳定 。 


表 2 隐 误 方法 和 显 式 方法 时 间 步 长 
Table BTiime step suited for explicit and 
\ implicit DG 


€ accuracy 


7.0x10-6 
5.0x10-4 


1.0x10-5 
1.0x10-3 


2.0x10-5 
2.0x10-3 


ime step(explicit) 


图 7 为 4 阶 精度 下 某 一 时 刻 的 涡 量 云图 ， 图 中 
涡 量 值 的 分 布 为 2; = 10x (i 一 1) 一 100, j= 1 …，21。 
图 8 和 图 9 分 别 为 从 同一 初始 流 场 开始 计算 ， 显 式 
方法 和 隐 式 方法 得 到 的 上 = 0.2 时 刻 壁面 摩擦 阻力 与 
压力 系数 系数 分 布 。 如 图 所 示 ， 隐 式 方法 的 计算 结 
果 与 显 式 方法 几乎 重合 ， 这 表明 了 隐 式 方法 的 正确 
性 . 图 10 和 图 11 分 别 表示 升力 系数 和 阻力 系数 随 
碗 代步 数 的 变化 图 , 当 非 定常 特性 保持 稳定 后 , 随 着 
涡 的 脱落 升力 系数 和 阻力 系数 均 呈 现 周 期 性 变化 规 
律 . 当 两 者 变化 保持 周期 性 时 , 提取 一 个 周期 的 峰值 
差 得 到 如 表 3 所 示 的 物理 参数 . 在 4 阶 精 度 (ps) 下 
与 Collist'7| 的 结果 进行 对 比 发 现 , 升力 系数 Cr 和 阻 
力 系数 Cp 的 峰值 误差 分 别 为 1.8% 和 6.0%, 表征 非 
定常 的 无 量 纲 参 数 斯 特 劳 哈 尔 数 5t 的 误差 为 1.1%， 
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这 表明 本 文 构造 的 非 定常 隐 式 DG 方法 的 精度 能 够 
达到 预期 值 . 


Level 1 2 3 4 5 6 7 8 9 10 11121314151617181920 
Oricity: 00 _90_80 -70-60-50-40-30-20-10 10 20 30 40 50 60 70 80 9010 


图 7 涡 量 云图 (ps) 


Fig. 7 Vortex contour (p3) 


0.30 


一 一 :6 一 :一 Explicit_Up 
0.25 一 一 全 一 一 Explicit Down 
Sd ed 


implicit_Up 
Implicit Down 


图 8 t= 0.2 时 刻 摩擦 阻 力 系 数 分 布 


Fig. 8 Skin friction coefficient distribution at t= 0.2 8 


1.5 © 


1.0| 优 
一 一 ' 避 .一 :一 Explicit up 
一 一 全 一 一 Explicit down 
0.5 一 一 会. 一 :一 Implicit up 
一 _A 


Implicit down 


3 


图 9 t= 0.2 时 刻 压力 系数 分 布 


Fig. 9 Pressure coefficient distribution at t= 0.2s 
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图 10 升力 系数 随 迭 代步 数 变化 图 


Fig. 10 Aift coefficient changing with iteration steps 
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表 3 不 同 精度 下 物理 参数 值 


Table 3 Physic parameter at different accuracy 


CD Cr St 
p1 1.3334 0.5991 0.1659 
p2 1.4432 0.7238 0.1672 
p3 1.4366 0.7380 0.1668 
CollisB29g 1.4107 0.6960 0.165 


结 


4 论 

本 文 针 对 非 定常 层 流 流动 的 隐 式 DG 方法 的 构 
造 进行 研究 , 主要 结论 如 下 ， 

1) 相 比 于 显 式 方法 ， 基 于 本 文 构造 的 非 定常 流 
动 的 隐 式 DG 方法 对 定常 流动 也 具有 较 高 的 求解 效 
率 , 且 实 际 精度 均 能 达到 设计 值 . 另外, 高 阶 壁 面 近 
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似 对 高 精度 CFD 方法 至 关 重 要 ， 

2) 非 定常 流动 的 隐 式 DG 方法 能 够 在 大 于 显 式 
方法 两 个 数量 级 的 时 间 步 长 上 保持 稳定 ， 表 明了 隐 
式 DG 方法 在 求解 非 定常 问题 时 具有 达到 较 高 求解 
效率 的 潜力 

3) 本 文 构造 的 非 定常 流动 的 隐 式 DG 方法 与 显 
式 方法 及 相关 文献 比较 均 有 较 高 的 吻合 度 ， 表 明了 
构造 方法 的 正确 性 ， 
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